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Abstract. The locations of massive stars 8Mq) within their host galaxies is reviewed. 
These range from distributed OB associations to dense star clusters within giant Hn regions. 
A comparison between massive stars and the environments of core-collapse supernovae and 
long duration Gamma Ray Bursts is made, both at low and high redshift. We also address the 
question of the upper stellar mass limit, since very massive stars (VMS, Mtnit ^> IOOMq) may 
produce exceptionally bright core-collapse supernovae or pair instability supernovae. 
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1.1. Star clusters 



C/3 ■ 1. Environments of Massive Stars 

1 Massive star formation in the Milky Way spans a broad spectrum, from dispersed, low 

Q-f intensity OB associations to concentrated, high intensity starbursts. Within a few hun- 
dred parsec of the Sun, high mass stars (Mj n i t ^ 8M©) are rather distributed, typically 
located in loose, spatially extended OB associations (de Zeeuw et al. 1999). A notable ex- 
ception is Orion OBI, which hosts the Orion Nebula Cluster (ONC), responsible for our 
i closest H II region. Further afield, large numbers of massive stars are associated with rela- 

tively intense bursts of star formation such as the high mass, compact clusters (Trumpler 
14, 16) within the Carina Nebula giant H II region. 

vo. 

It is generally accepted that the majority of stars form within star clusters (Lada & Lada 
2003), although recent evidence suggests star formation occurs in a continuum of stellar 
ON ' densities (e.g. Evans et al. 2009). Nevertheless, given their short-lifetimes (3-50 Myr) 
only a few percent of massive stars appear genuinely 'isolated' (de Wit et al. 2005) such 
that they either tend to be associated with their natal cluster or are plausible runaways 
. . . from ijj]. 

According to Weidner & Kroupa (2006), there is a tight relation between cluster mass, 
and the most massive star formed within the cluster, although this remains controversial 
(Calzetti et al. 2010, Eldridge 2012). Examples of well known star clusters spanning a 
range of masses are shown in Table Q] all of which are sufficiently young (< 1 — 2 Myr) 
that the most massive stars have yet to end their lives. We include the most massive star 
in each cluster, which increases towards the highest mass clusters. 

If there is a relation between a cluster and its most massive star, the galaxy-wide 
stellar initial mass function (IMF) will also depend upon the cluster mass function and 

t Runaways may be ejected from their cluster either dynamically during the formation process 
or at a later stage after receiving a kick following a supernova explosion in a close binary system. 
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Table 1. Selected young star clusters spanning a range of masses, M c \, for which (initial) 
masses of the highest mass stars, M*,i n u, have been determined. 



Cluster 


M cl /M 


Ref 


Star 


Af», init /M Q 


Ref 


p Oph 


~ 10 2 


a 


p Oph Source 1 


9 


a 


ONC 


1.8 x 10 3 


b 


6 1 Ori C 


39±6 


c 


NGC 3603 (HD 97950) 


~ 10 4 


d 


NGC 3603-B 


166±20 


e 


R136 (HD 38268) 


5 x 10 4 


c 


R136al 


320±™ 


c 



(a) Wilking et al. (1989); (b) Hillenbrand & Hartmann (1998); (c) Simon-Di'az et al. (2006) 
(d) Harayama ct al. (2008); (c) Crowthcr et al. (2010) 

the range of cluster masses, as set out by Pnamm-Altenburg et al. (2007). In normal star- 
forming galaxies the cluster mass distribution follows a power law with index —2, albeit 
this is truncated at high mass depending upon the rate of star formation (Gieles 2009). 
Consequently, similar absolute numbers of stars are formed in low mass {M c \ ~ 1O 2 M0), 
intermediate mass (~ 10 3 M©) and high mass (~ 10 4 M Q ) clusters, while high mass stars 
should be rare in the former. This is not always the case, since star formation in some 
nearby dwarf irregular starbursts is strongly biased towards a few very high mass clusters 
(e.g. NGC 1569, Hunter et al. 2000). 

1.2. Hu regions and star formation rates 

Because of the (universal?) Salpeter IMF slope, the overall statistics of massive stars in 
galaxies will be heavily biased towards 8-20 Mq (early B-type) stars. However, the most 
frequently used indicator of active star formation is nebular hydrogen emission (e.g. Ha) 
from gas associated with young, massive stars. The Lyman continuum ionizing output 
from hot, young stars is a very sensitive function of temperature (stellar mass), such that 
one 03 dwarf (~75 Mq) will emit more ionizing photons than 25,000 B2 dwarfs (~9 Mq, 
Conti et al. 2008). Therefore, Hn regions are biased towards high mass (O-type) stars 
with >20 Mq since B stars will produce extremely faint Hn regions. 

Beyond several Mpc, current sensisitivies limit detections of Hn regions to relatively 
bright examples, involving several ionizing early O-type stars (Pnamm-Altenburg et al. 
2007) . Still, the Ha luminosity of bright H n regions can be converted into the correspond- 
ing number of Lyman continuum ionizing photons, for which the number of equivalent 
07 dwarf stars, N(07V), serves as a useful reference (Vacca & Conti 1992), as indicated 
in Table[2j Kennicutt et al. (1989) have studied the behaviour of the Hn region luminos- 
ity function in nearby spirals and irregular galaxies. Early-type (Sa-Sb) spirals possess 
a steep luminosity function, with the bulk of massive star formation occuring in small 
regions ionized by one of a few O stars, plus a low cut-off to the luminosity function. 
Late-type spirals and irregulars possess a shallower luminosity function, in which most 
of the massive stars form within large Hn regions/OB complexes, for which 30 Doradus 
in the LMC serves as a useful template. For example, although the LMC contains con- 
siderably fewer Hn regions than M31 (SAb), it contains ten Hn regions more luminous 
than any counterpart in M31 (Kennicutt et al. 1989). 

The integrated nebular Ha luminosity of a galaxy is widely used as a proxy for the rate 
of (near-instantaneous) star formation (Kennicutt 1998), although conversions into total 
star formation rates (SFR) rely upon the adopted stellar mass function and evolutionary 
models for single and binary stars (e.g. Leitherer 2008). In addition, since the youngest 
star forming regions are deeply embedded, the combination of gas (Ha) and dust (24/im 
continuum) provide a more complete SFR indicator (Calzetti et al. 2007), although the 
situation is more complicated for galaxies with low SFR (e.g. Pnamm-Altenburg et al. 
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Table 2. Examples of nearby Hn regions, spanning a range of luminosities (adapted from 
Kennicutt 1984), for an assumed 07V Lyman continuum ionizing flux of 10 49 ph/s). 



Region 


type 


galaxy 


Distance 


Diameter 




IN I I V ) 








(kpcj 


\P C ) 


(ergs 1 




Orion (M42) 


Classical 


Milky Way 


0.5 


5 


lxlO 3 


<l 


Rosette (NGC 2244) 


Classical 


Milky Way 


1.5 


50 


9xl0 3r 


7 


N66 


Giant 


SMC 


60 


220: 


6xl0 38 


50 


Carina (NGC 3372) 


Giant 


Milky Way 


2.3 


300: 


l.5xl0 39 


120 


NGC 604 


Giant 


M33 


800 


400 


4.5xl0 39 


320 


30 Doradus 


(Super)giant 


LMC 


50 


370 


l.5xl0 40 


llOO 


NGC 5461 


(Super)giant 


M101 


6400 


1000: 


7xl0 40 


5000 



2007). In addition, Ha-derived star formation rates differ from FUV continuum diagnos- 
tics for dwarf galaxies (Lee et al. 2009b), while FUV indicators closely match the local 
ccSNe rate (Botticella et al. 2012). 

1.3. 30 Doradus: Template extragalactic giant Hll region 

30 Doradus, the brightest star forming complex within the Local Group, provides a useful 
template for extragalactic 'supergiant' Hll regions (Kennicutt et al. 1995, Table[2]). The 
30 Dor nebula is shown in Fig.[T]and spans an angular size of ~ 15' x 15', corresponding 
to a linear scale of 220 x 220 pc at the distance of the LMC. Consequently, individ- 
ual stars may be studied in detail (e.g. Evans et al. 2011). Walborn & Blades (1997) 
identified five distinct spatial structures within 30 Dor, (i) the central 1-2 Myr cluster 
R136; (ii) a surrounding triggered generation embedded in dense knots (< 1 Myr); (iii) a 
OB supergiants spread throughout the region (4-6 Myr); (iv) an OB association to the 
southeast surrounding R143 (^5 Myr); (v) an older cluster containing red supergiants to 
the northwest (10-20 Myr). 30 Dor would only subtend 1.5" at a distance of 30 Mpc, so 
care should be taken for nebular-derived ages of stars within extra-galactic H II regions 
(e.g. Leloudas et al. 2011). 

2. Environments of supernovae and gamma-ray bursts 

2.1. Hll regions and core-collapse SNe 

Turning to studies of the environments of supernovae, locally neither type II nor type 
Ib/c supernovae are associated with ongoing star formation. Specifically, Smartt (2009) 
examined the host environment of a volume limited {cz < 2,000 km/s), statistically 
complete sample of ccSNe, of which from 20 type II SN were located in bright H n 
regions. A number of type II SN were located in loose associations, with two in older 
clusters (e.g. SN2004am, II- P, in M82), while only 1 of 10 type Ib/c SN from Smartt 
(2009) was in a large star forming region (SN2007gr, Ic, in NGC 1058), albeit spatially 
offset from regions of H II emission. 

Anderson & James (2008) took a different approach, studying the assocation between 
ccSNe and H II regions within (mostly) bright spirals, whose recession velocities extended 
up to cz = 10,000 km/s. In common with Smartt (2009), Anderson & James (2008) did 
not find type II SNe associated with H n regions, concluding that the "type II progen- 
itor population does not trace the underlying star formation". In contrast, Anderson & 
James (2008) found that type lb, and especially Ic ccSNe were spatially coincident with 
(presumably bright) Hll regions. 

Let us consider the typical duration of the H II phase in young, isolated clusters. Wal- 



12 



Paul A. Crowther 




30 Dor 

(LMC) 



Carina 
Nebula 



N66 (SMC) 




150 pc 



Figure 1. Three Local Group giant H n regions shown on the same physical scale: Carina Nebula 
(Milky Way, ESO/WFI, 60' x 30'), 30 Doradus (LMC, ESO/TRAPPIST, 20' x 20'). N66 (SMC, 
HST/ACS, 4.7' x 4.7'). 30 Doradus hosts multiple stellar generations (Walborn & Blades 1997) 
but would only subtend 1.5" at a distance of 30 Mpc, so care should be taken for characteristic 
ages of stars within extra-galactic H n regions. 



born (2010) compared the properties of young star clusters, revealing an association with 
a Hn region only for the first ~2-3 Myr, after which the gas has been dispersed (e.g. 
Westerlund 1, Clark et al. 2005). Therefore, one would not expect ccSNe to be spatially 
coincident with isolated H II regions unless the mass of the progenitor was sufficiently 
short for its lifetime to be comparable to the gas dispersion timescale. This is illustrated 
in Figure[2^a) where we compare the lifetime of the most massive stars in clusters (masses 
according to Eqn. 10 from Pflamm-Altenberg et al. 2007), adopting stellar lifetimes from 
Ekstrom et al. (2012), with an estimate of the duration of isolated Hn regions (adapted 
from Walborn 2010). This naturally explains the lack of any association between type II 
ccSNe and Hii regions for both Smartt (2009) and Anderson & James (2008). 

How, then, can one explain the empirical association between type Ib/c SNe and Hn 
regions in the Anderson & James (2008) study? These either arise from very massive stars, 
which would be inconsistent with Smartt (2009), or more likely we have to appreciate 
that not all massive star formation occurs within isolated, compact star clusters. 

Late-type spirals and irregulars, which form the majority of Anderson & James' host 
galaxy sample, host large star forming complexes, up to several hundred parsec in size, 
involving (super)giant Hn regions. These are ionized by successive generations of star 
clusters, separated by a few Myr (Tabled, with a total duty cycle of ^10 Myr. There- 
fore, a massive star exploding within such an environment as a SN after 5-10+ Myr 
would still be associated with a bright Hn region, as illustrated in Fig. HJb), even if its 
natal star cluster had cleared the gas from its immediate vicinity. Resolving the location 
of the ccSNe within the region would be especially difficult at larger distances. Recall 
that the average distance of galaxies within the Anderson & James (2008) sample was 
~32 Mpc, and that their study was based upon moderate resolution ground-based Ha 
imaging. Typically higher spatial resolution datasets were employed by Smartt (2009), 
which together with a lower host distance (^27 Mpc maximum for an adopted Hq = 75 
km/s/Mpc) enabled a higher spatial inspection of the SN environment (recall Fig. []}. 
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Figure 2. (a) Schematic comparing the lifetime of the most massive star in a cluster (accord- 
ing to Pfiamm-Eltenberg et al. 2007) and isolated Hn regions (adapted from Walborn 2010). 
Core-collapse SNe should only be associated with isolated Hn regions for very massive progeni- 
tors; (b) as (a) except for (super)giant H n regions, whose 10-20 Myr lifetimes should imply an 
association with ccSNe except for only relatively low mass, long lived (type II-P) progenitors. 



2.2. Wolf- Ray et stars and ccSNe 

Main-sequence O stars may preceed ccSNe by up to 3-10 Myr, whereas Wolf-Rayet stars, 
their evolved descendents, should preceed (type Ib/c) SNe by a timescale that is an order 
of magnitude shorter (Crowther 2007). Therefore, comparisons between the environment 
of Wolf-Rayet stars and ccSNe provide information upon whether the former are plausibly 
the parent population of the latter (e.g. Leloudas et al. 2010), since lower mass close 
binaries might dominate type Ib/c SNe statistics (Smith et al. 2011). In addition, we can 
also compare the environment of Wolf-Rayet stars in their host galaxies (e.g. associated 
H II regions) with those of ccSNe. 

Since the Wolf-Rayet content of the Milky Way is highly incomplete due to foreground 
interstellar dust, let us instead turn to local external galaxies. From the LMC Wolf-Rayet 
catalogue of Breysacher et al. (1999), 58% (78/134) lie within OB associations, while 
83% (112/134) lie within catalogued Hn regions. Of course, the LMC is not particularly 
representative of the local star forming galaxy population, since star formation is largely 
confined to several giant H n regions (Kennicutt et al. 1995) and faint H n regions would 
not necessarily be identified in distant galaxies. Neugent & Massey (2011) have re-assessed 
the Wolf-Rayet content of M33 (Scd), revealing a total of 206 stars which they argue 
is complete to ^5%. As for the LMC, the majority of WR stars (80%) reside in OB 
associations. Further afield, a comparison between the recent Wolf-Rayet photometric 
survey of the Scd spiral NGC 5068 (Bibby & Crowther 2012) with Ha images reveals 
that 50% of the Wolf-Rayet candidates lie in bright or giant Hn regions, while 25% are 
associated with faint H n regions and 25% lie away from any nebulosity. 

Recalling Sect. 12. ll type Ib/c SN are rarely associated with Hn regions. This suggests 
that most do not result from young, massive Wolf-Rayet stars, arising instead from lower 
mass close binaries in which the ccSNe arises from the H-deficient, mass-losing primary 
(Fryer et al. 2007). In contrast, the preference of type Ic SNe for Hn regions suggest 
that massive Wolf-Rayet stars are realistic progenitors in such instances. Close binaries 
in such a scenario ought to mimic the masses, and in turn, the lifetimes of type II ccSNe 
progenitors, which positively shy away from Hn regions. 
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Table 3. Summary of expected association between Hn regions and ccSNe/long GRBs in 
different host galaxies (following Kennicutt et al. 1989, Gieles 2009). 



Host 


SFR 


Cluster range 


Characteristic 


SN-Hn 


Example 






(M Q ) 


Hn region 


association? 




Spiral (Sab) 


Low 


10 a-4 


Isolated 


No (all types) 


M31 


Spiral (Scd) 


High 


10 a-6 


Giant 


Yes (Ib/c), No (II-P) 


M101 


Irr 


Low 


10 a-4 


Isolated 


No (all types) 


SMC 


Irr 


High 


10 2- 6 


Giant 


Yes (Ib/c, GRB), No (II-P) 


NGC 1569 



2.3. ccSNe, long GRBs and host galaxy types 

Mindful of the spatial resolution issue, let us turn to the Kelly et al. (2008) study of SNe 
locations with respect to the continuum g'-band light from their low redshift (z < 0.06) 
host galaxies. Kelly et al. revealed that Ic SNe are much more likely to be found in the 
brightest regions of their hosts than lb or II SNe. An earlier analysis of high redshift 
galaxies by Fruchter et al. (2006) revealed that long GRBs (< z > = 1.25) were also 
strongly biased towards the brightest pixel of their hosts, in contrast to core-collapse 
SNe (< z > — 0.63, presumably mostly type II-P) which merely traced the light from 
their hosts. Kelly et al. (2008) concluded that if the brightest locations correspond to 
the largest star-forming regions, type Ic SNe (and long GRBs) are restricted to the most 
massive stars, while type lb and especially type II-P SNe are drawn from stars with more 
moderate masses, results in common with Anderson & James (2008). 

However, one significant difference between the low-redshift SN Ha study of Anderson 
& James (2008) and the high-redshift GRB study of Fruchter et al. (2006) is that hosts 
of the former are relatively high mass, metal-rich spirals, while those of the latter are 
low mass, metal-poor dwarfs. In normal disk galaxies the number of stars forming across 
the mass distribution of star clusters is relatively flat, albeit with a cut-off linked to 
the star formation intensity (Gieles 2009). The star cluster mass function is repeated in 
nearby dwarf galaxies (Cook et al. 2012), but galaxy- wide triggers may induce intense, 
concentrated bursts of star formation, leading to disproportionately numerous massive 
star clusters (Billett et al. 2002 jj We have attempted to set out the potential association 
between H n regions, ccSNe and long GRBs in Table [3] for star forming spirals and 
irregulars, based upon the above arguments, although exceptions are anticipated (and 
subject to uncertainties regarding the main progenitors of type Ib/c SNe). 

Relatively massive, metal-rich galaxies would represent the primary site of all star 
formation for the sample of Fruchter et al. (2006), resulting in (type II-P) ccSNe unasso- 
ciated with the brightest regions in their hosts. Yet, when localised starburst activity does 
occur, it is very intense (Billett et al. 2002), leading to very massive clusters, and in turn 
large numbers of high mass, metal-poor stars, a subset of which would be progenitors of 
the long GRBs witnessed by Fruchter et al. (2006). 

3. Upper Mass Limit 

The lower limit to the mass of stars is relatively well known (e.g. Burrows et al. 1993), 
yet establishing whether there is a corresponding upper mass limit has proved elusive 
(Massey 2011). In part, this is because obtaining robust masses for VMS is extremely 
challenging, and in part because of the scarcity of star clusters that are sufficiently 

f Of course, not all dwarf galaxies are starbursting. Within the local volume (<11 Mpc) only 
a quarter of the star formation from dwarf galaxies is formed during starbursts (Lee et al. 2009a) 
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nearby, young and massive for their most massive stars to be studied in detail. Up until 
recently, a mass limit of <~ 15OM has been commonly adopted, based upon a near-IR 
photometric study of the Arches cluster (Figer 2005). However, it is well known that 
the temperature of hot, massive stars is rather insensitive to optical/IR photometry. 
Spectroscopic analysis is required for robust temperatures and in turn luminosities, from 
which stellar masses are derived. 

3.1. R136 stars 

The situation is especially difficult for the brightest main-sequence members of the most 
massive young clusters, which possess unusual (emission line) spectral morphologies, rem- 
iniscent of Wolf-Rayet stars (e.g. Drissen et al. 1995). The mass-luminosity relationship 
for main-sequence VMS is relatively fiat, L oc M 1 - 5 (e.g. Crowther et al. 2012), so in- 
ferred masses are particularly sensitive to temperature, M oc T c ^ . Recent advances in 
atmospheric models for stars with dense stellar winds has led to an upward revision to 
the temperatures of such stars, by ~25%, corresponding to as much as an 80% increase 
in the resulting mass. Fortunately, several very massive, double-lined eclipsing binaries 
have been identified within the past few years, including the Wolf-Rayet binary NGC 
3603 Al (Schnurr et al. 2008), permitting an independent check on spectroscopic results 
for similar systems. 

R136, the central ionizing cluster of 30 Dor, has both a very high stellar mass (^55,000 
Mq) and a sufficiently young age (1-2 Myr) for its most massive stars not to have under- 
gone core-collapse. Previous estimates of their stellar masses, based on conventional O 
star calibrations, implied 120 - 155 M Q (Massey & Hunter 1998). Schnurr et al. (2009) 
searched for close binaries among the visually brightest members, but none revealed ra- 
dial velocities, with the possible exception of R136c. Still, their near-IR integral field 
datasets provided spatially resolved spectroscopy of individual stars within R136, which, 
together with archival UV/optical spectroscopy and AO-assisted photometry permitted 
a reassessment of their stellar masses. Spectroscopic analyses together with new evolu- 
tionary models for VMS enabled Crowther et al. (2010) to revise their (current) stellar 
masses upward to 135-265 M . Initial masses of 165-320 M Q were inferred, adopt- 
ing standard main-sequence mass-loss rates for VMS (Vink et al. 2001) which closely 
matched spectroscopically-derived values, and were reinforced by the close agreement 
between spectroscopic and dynamical masses obtained for NGC 3603-Al. Overall, R136 
supports the trend that higher (initial) mass stars reside within the most massive star 
clusters set out by Weidner & Kroupa (2006). However, statistics of high mass clusters 
for which accurate stellar masses have been determined remain very poor. 

3.2. Pair instability supernovae 

Based upon their re-assessment of the most massive stars in R136 and other young, high 
mass clusters (Arches, NGC 3603), Crowther et al. (2010) concluded that their stellar 
content was consistent with a revised upper mass limit of ~ 300M Q . Regardless of the 
physical origin of this limit, such high initial masses raise the prospect of extremely lumi- 
nous core-collapse SNe (Waldman 2008) or even pair-instability SNe (Heger & Woosley 
2002). Models have recently been calculated for the post-main sequence evolution of VMS 
spanning a range of metallicities (N. Yusof, these proc). From these, would appear that 
the VMS in R136 will end their lives as core-collapse SNe, with lower metallicity (SMC- 
likc) required to reduce mass-loss rates sufficiently for pair-instability SNe, as has been 
proposed for SN 2007bi (Gal- Yam et al. 2009). However, details remain very sensitive to 
mass-loss prescriptions for the post-main sequence evolution (e.g. Crowther et al. 2012). 
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